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- - - Abstract 

cn ; 

. Quadratic Lagrangians are introduced adding surface terms to a free 

I particle Lagrangian. Geodesic equations are used in the context of the 

Hamilton-Jacobi formulation of constrained sysytem. Manifold structure 
induced by the quadratic Lagrangian is investigated. 
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1 Introduction 



J> ■ There are mainly two basic methods to investigate constrained systems. First 

■ one is the method initiated by Dirac [1]. This method is based upon the clas- 

O ■ sification of constraints and consistency conditions derived from the the vari- 

^ ! ations of the constraints. The symplectic structure is established by defining 

I Dirac brackets. The second method is the treatment of the constrained sys- 

O ■ tems by Caratheodory's equivalent Lagrangian method [2], [3], [4]. Since this 

method leads us to the Hamilton-Jacobi equations it will be called Hamilton- 
I . Jacobi (H J) approach. Equivalent Lagrangian method [2] leads us to a set of 

{^j;)! Hamilton-Jacobi equations [3], [4]. Despite of many attempts to clarify the many 

aspects of Hamilton-Jacobi formalism (HJ) we still have some subtle ones. The 
' main difference between Dirac's formulation [1] and (HJ) approach [5], [6], [7] 



' is the fact that the second one is multi Hamiltonian approach depending on the 

rank of the Hessian matrix. On the other hand the boundary conditions [8], [9], 
[10], [11] play an important role for constrained systems [12] and they are cru- 
cial for (HJ) formulation [2], [13]. The presence of many Hamiltonians, which 
are constraints having a special forms, suggests us to investigate the possibil- 
ity to describe the superintegrable systems [14] and to obtain some integrable 
geometries. 

In this paper we start with a free Lagrangian theory and we add a surface term 
involving the components of the angular momentum. By a suitable choosing of a 
Lagrange multipliers we obtain a quadratic Lagrangian which can be singular or 
not. The aim of this paper is to study quadratic, singular Lagrangians, using the 
(HJ) formulation. This paper is arranged as follows: In Sec. 2 (HJ) formalism 
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is briefly explained. In Sec. 3 the method of obtaining quadratic Lagrangians 
is presented and two non-singular quadratic systems are studied. In Sec. 4 the 
singular case is analyzed using the (HJ) formalism. Conclusions are presented in 
Sec. 5. In Annex some geometrical properties of the obtained metrics are shown. 



2 Hamilton-Jacobi formalism 

Let us assume that the Lagrangian L is singular and the Hessian matrix r-, 
a,&=l,2,..,n — r, has rank n-r.The canonical Hamiltonian Hq is defined as 

The other Hamiltonians are obtained from the definitions = ^ |ga=wa ^^^1 
they have the forms 

where a, (3 — n — r + 1, ..,n; a — 1,2, .., n — r. Thus, one arrives at the following 
set of Hamilton-Jacobi partial differential equations 

^ + ^"(^^9"'^>'^ = ^)=0 (2) 

The following total differential equations lead us to determine 
the Hamilton-Jacobi function S{if^,q°'): 

r)H' 

dq^ = ^dr (3) 

dpa 

dpa = -^dr (4) 
r)H' 

dp a = -^dr (5) 
Solutions of these equations S{q^ as 

dz^{-H^+p,—^)dr, (6) 

dpa 

where z — S{t°', g"). Since the equations of motion are total differential equations, 
one should check the differentiability conditions, consistency conditions in Dirac 
formulation. It is proved that the system of equations (3) is differentiable if and 
only if the following Poisson brackets equations are satisfied : 

[<,^^] = 0. (7) 

In general some of these Poisson brackets may not vanish identically. In such 
cases one should include the new functions as "new" constraints. Poisson brackets 
of these constraints with should be considered also. This procedure should 
continue until all Poisson brackets vanish identically. This approach was studied 
in detail in references [3], [4]. 
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3 The method 



Let us assume that a given Lagrangian L{(f,q^) = \aij(fq^ admits a set of inde- 
pendent constants of motion denoted by Lj, i = 1, • • • n. Our aim is to treat the 
new Lagrangian L — L -\- A*Lj as a quadratic Lagrangian expresses as 

^' = \gijq'q'- (8) 

Let us assume that the matrix gij is non-singular. 

Treating gij as the " metric tensor" of a manifold the geodesic equations 

^ 2^ ^ dq^ ^ ag- dq^ ^ 

are equivalent to the Euler-Lagrange equations of (8). Here g"-^ is the inverse 
of gab- In this formulation the canonical Hamiltonian of the system is expressed 
as H — ^g"'^PaPb- In this way the motion of a particle is described by geodesic 
equations of a Riemannian manifold with metric tensor g^j. The above motion is 
treated inside of (HJ) formalism and the geometrical properties of the obtained 
metrics g^j will be investigated. In the following we will add to a free Lagrangian 
the components of the angular momentum. 

l.As the first example of a non-singular quadratic system, let's consider the 
following Lagrangian 

L^^{x^ + f) + Xs{xy-yx), (10) 



which can be expressed ci\ where 

1 

(11) 



ttij — 



/I -y \ 

Ola; 
\-y X / 



This Lagrangian corresponds to the motion of a particle moving on a plane such 

that Lz is constant. 

Hamiltonian has the following form 

He = ^%pPi,Pp where //, p = 1, .., 3 (12) 

and a~p is the inverse of a^p. Prom Hamiltonian we find 

dx = aipPpdt, dy = a^pPpdt, d\^ = a^pdt 
dPx^^{^)p^Ppdt, dpy = f{^)p^ppdt, dpx,^0 ^^^^ 
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In order to find equation of motion we must solve (13). However we will prefer 
geodesic equations that produce easier equations and equivalent to the (13). The 
geodesic equations corresponding to (11) are 



X = 



y 
As 



'-■^[xy - yx), 

+y^(xy-yx), 
^{xx + yy). 



(14) 



The aim is to solve the system of equations given by (14). Notice that (14) implies 
that 

A3 ^ -2{xx + yy) 
A3" 



which admits solutions as 



A. 



/ 



C 



+ y"^ 

where C is constant. Taking into account that 



x 

y 



X 

y 



and denoting 



dt, 



X 
X 



(15) 
(16) 



= k 



y 



we divide our problem into two parts. The first one corresponds to the case when 
A; < 0. We obtain the solutions as 

x{t) — ci sm{Vkt) + C2 cos{Vkt), y{t) — C3 sin{Vkt) + C4 cos{Vkt). (17) 

If we impose some initial conditions one can find specific solutions of equations. Let 
us assume that ci — 0, — 1, C3 — 1 and C4 = 0. Prom (17) we find that 
x{t) — cos(a;t) and y{t) = sin(a;t) where uj = \fk. Also from equation (16) we 
obtain that A3 = Ct + i? where B is a constant. On the other hand Vy ^iid 
PA3 are given as follows 



(18) 




1 





-»\ 







1 


X 




-y 


X 


) 





PA3 



If we solve (18) we find 

'Px = sm{ujt){—uj — C) Py = cos{ujt){uj + C) 
In the case of A; > the solutions are given by 



(19) 



(20) 



where q, i = 1, .., 4 are constants. As in the previous case A3 is given by (16). Using 
(20) and imposing A = C3C2 — C1C4 7^ we obtain the following equation (c| + 
c|)x^ + (c? + c|)t/2 - 2xy{c2C4 + C1C3) = A^. 
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The metric from (11) is conformaly flat having its Ricci scalar R — ■^2^- 
2. Let us add now two components of the angular momentum to the La- 
grangian of a free three-dimensional particle. We obtain 



L — -(x^ + 2/^ + z'^) + Ai(|/i; — zy) + X2{zx — xz), 



(21) 



and from (21) we identify the metric 



aij — 



/ 1 

1 



-z 

\ z 





1 

y 



-X 





—z 

y 






z \ 


-X 



J 



(22) 



X 



In this case the geodesic equations have the following forms 

■[2\i{yz — zy) + 2X2{zx — xz)] 

y = -jT—^—^i'^My^ - ^y) + '^>^2{zx - xz)\ 

^[2\i{yz - zy) + 2X2{zx - xz)] 



z — 



X 


_|_ y2 




y 




x^ + y^ 




z 





2-2 _|_ ^2 _|_ 



Ai 

A2 



^2 _|_ ^2 _|_ ^2 



[-Ai 



\x\\ 



,x?z + z^z 



+ yy) + A2yi 



zx — xz 



)] 



,zy - yz 



z 



^) - A2( 



y^z 



z^z 



XX] 



(23) 
(24) 
(25) 

(26) 

(27) 



2;2 _|_ ^2 _|_ ^2 

These equations form a system of nonlinear ordinary differential equation and to 
solve this system we should notice that 2 = ^ = 2. Assuming that the above 

ratio is a negative constant k we obtain 

x{t) = ci s\n{\/kt) + C2 cos(v^t) (28) 

y{t) = C3 sin(^t) + C4 cos{Vkt) (29) 
-2(^) = C5 sin(v^t) + cg cos(V^i), (30) 
where —k' = k and the constants q, i = 1...4 are subjected to the restriction 



that — ^ ~2^^2^^2^^~ is a constant. Note that if /c = 0, both x and y must be 0, 
then X and y are both linear function of t. Ai and A2 can be easily obtained from 
(26) and (27). For the general case the expressions involved are very complicated. 
Instead of doing this we will present a particular solution. If we assume that 
C2 — ^ , cs — 1, — ^ and others are , we obtain 



X 



(t) — -—cos{Vkt), y{t) — sin{Vkt) and z{t) — -—cos{Vkt). (31) 
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Using (31), (26) and (27) we obtain 

Xi{t) = (32) 
K{t) = v^Ai(t) \ ; A2(t). (33) 

COS( V J 

The solution of (32) is 

Ai(t) = ht + h, 

Hi) = jhln{l + tan^{y/kt))V2^ + \hln{cos{^t))V2^ (34) 
+ ^kisin{2^/kt) + (^ - ^ cos(2v^t))A;3 + ^2 + \ht 

where ki,i — 1, ■ ■ ■ A are constants. As before Pi have the forms Pi — aijqL 



4 The singular case 

As an example of singular quadratic system, let's consider the Lagrangian given 
by ^ 

L — -(x^ + + z^) + Xiiyz — zy) + A2(2;i; — xz) + X'i^xy — yx) (35) 

which is expressed in compact form as L = ^aijq'^q^, where the matrix a^j has the 
form 





1 











z 


-y 









1 





—z 





X 












1 


y 


—X 












—z 


y 















z 





— X 













I 


-y 


X 














J 



(36) 



Since the rank of the matrix is 5, this is a singular system. By using Pi — ^ 
we find 

Px + yXs ^ X + zX2, py- xXs^ y - zXi,p^^ a3iq\ px^^aaq', p\^ ^ anq\{'il) 
From (37) we find 

■ _ Py zy^ +Pz z'^y+Px j^z^+xyzp^- 3:ypx^ -x'^PX^ 

. _ Pzz'^+Pyyz+px-^^y-xX2-xzpx /„„n 
^ ~ zix^+y^+z'^) ' ^ y'^^) 

_ -pyZ^+z^x\3+pzZ^y-pxj^z^+xyzpx+x"'zX[i-zx^Py+zx\3y'^~xypx2-x^PXi 
~ zix^+y^+f) ^ ' 

^ _ -px2Z^-P\2y^z^Px+zy^px+z^y\3+y^zX3-pzXZ^-pxiXy-pyzyx+x^yz\3 
— zix^+y'^+z^) 
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Using (38), the canonical Hamiltonian has the following form 



where h,,} is the inverse of 



yp\2 + xpx^ 



) where fj,,p—l,..,5 (39) 



/ 1 











z 


\ 





1 





—z 













1 


y 


—X 







—z 


y 










V ^ 





—X 








/ 



In (HJ) formalism the Hamiltonians to start with are 

Hq = Po + ^b^pPf^Pp, = PX3 + 

The total differential equation corresponding to (41)are 



dx 


= bipPpdt, 


dXi 


= b^pPpdt + 


dpxi 


= 0, dpx^ -- 


dpy 


2^ dy ' 


dpz 


-\db-l 
2^ dz ' 



-( 

z 



jpdt — ^^dXs, 
jPuPpdt + ( 5 jdXs 



P\i 



(40) 



(41) 

(42) 
(43) 

(44) 

(45) 

(46) 



The next step is to find the variations of Hq and H^. Using the fact that b~p is a 
symmetric matrix and taking into account the previous equations we obtain after 
some calculations that both variations are zero. We conclude that the system 
is integrablc and we observe that = zrxLi which is identically zero. The 
geometrical properties of the 5X5 metrics are presented in Annex. 



4.1 Curved space generalization 

Let us start with a free three dimensional Lagrangian 

L= i(x2 + ?/2 + z2) 

If we impose z"^ = 1 — x^ — and take its time derivative, we obtain 
z = So, (47) becomes 



(47) 



(48) 
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where 



9ab — ^ab + 



u 



(49) 



where z'^ = u, a,b = 1,2. The above metric admits three invariants [15] 

Lx = -Vupy, Ly = ^/upa:, = xpy - yp^, (50) 

Adding those three constants of motion to (48 ) we obtain the extended La- 
grangian in a compact form as L = k^'ijQ^Q^, where a^j is has the following form 



/ 1 + ^ 



xy 
u 



y 



X 



u 



\ 



-y 



xy _Ey. a: i Z:^ 

■^-Ju 

xy 













X 





/ 



(51) 



The determinant a^j is so is singular symmetric matrix and the rank of it is 
4. In order to find the canonical Hamiltonian we need x, y, Ai and A2 in terms 
of Px-,Py, Pai, P\2 and A3 .Using 

dL 

Pi = — = aijq^ where i,j = 1, .., 5, 



we obtain 



X — 



(52) 



l—x-^—y' 



A 



■ _ (pa 1 -PA 1 A2 ) \/i -^^ -y^ 

" 1— i^— y2 

_ x^\3+x^ypx -x^py +x^pyy'^ +xy^ A3 -2: A3 \-xy^ X-x^ -y'^+xy^px 
Ai — / ^ 

^yPi^-PAi\/l"-a:^-i/^i/^-2pyj/^+PAi\/l-a:^^-y^+;/'*Pi/+Py 

_ PA2 V 1- -y^ (x^ -1) -p:r +2a:^pi, -yA3 +a:^yA3 - J^Px +y^Pj! -y^x^Px 
^ ~ ^Jl-x^ -y^{l-x^-y^) 

y^ A'a -xy^py +xypx^ y/l-x^-y'^ -x^yPy+xypy 

Vl-x2-2/2(l-a;2-y2) 



(53) 



SO, the canonical Hamiltonian is given by 



He = -^gJPuPp + A3(pa3 + 



ypA2 + ^pai 

— x^ — y^ 



) where /i, p — 1, .., 4 



(54) 



where 



/ 1 + 



xy xy_ x2_ 

-i-v^ 











(55) 
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Since the matrix g^p is singular, we obtain a constrained system. In (HJ) formal- 
ism we have two Hamiltonians 

Ho^Po + ^g^^pPi^Pp, = pxs + ^y==^=^ (56) 



The total differential equations of (56) are given by 



dx = QipPpdt, dy = g2pPpdt, dXi = g-^pdt + -^j= 
d\2 = gXpdt + --j=J^=^dXz, 

The variations of and are 0. Geodesic equations corresponding to are 
given by 

yu — y{y^ — l)x^ — 2xy'^xy + y{x'^ — l)y^ (58) 

xu — x{y'^ — l)x^ — 2yx^xy + x{x^ — l)y^ (59) 

Aim = -2x{y'^ -l)x\i-2y{y'^ -l)x\2 + 2x^yy\i + 2xy^y\2 (60) 

A2M = -2x(a;^ - l)yAi - 2y{x^ - l)y\2 + 2x^yx\i + 2x?/^i;A2 (61) 

where u—l — x^ — y^. After some calculations we obtain a solution 
as 

x{t) — cos{kt),y{t) — sin{kt), 
Ai(t) = ie-^+'=*cos(H)C4 - |cos(H)C3e(-i/2— fc*) + Ci, (62) 

A2(t) = (C2k - C3e-^-'=*Sin(H) + g4e(-^+''*)sin(fct) _ 

The geometrical properties of the corresponding metrics are presented in An- 
nex. 



5 Conclusions 

In this paper we illustrated the importance of the surface terms in finding some 
integrable geometries in three, four and five dimensions. The central role was 
played by the components of the angular momentum. We added a surface term 
A3L2 to a free two-dimensional Lagrangian and we obtain a three dimensional 
geometry which is conformal flat and has a non-vanishing scalar curvature. Wc 
solved the geodesies equations and we found non-trivial solutions. We repeated 
the procedure and we investigated the case when two components of the angular 
momentum were added to a free three-dimensional Lagrangian.In this manner 
we obtain a five-dimensional metric admitting three Killing vectors. A singular 
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system was obtained if all components of the angular momentum were added. We 
treated this system within (HJ) formalism and we obtained three types of four 
dimensional metrics admitting three isometrics. The procedure was extended to 
the curved space and the two-dimensional sphere was investigated in details. 
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7 Annex 

In the following the Christoffel components, Ricci scalar and the Weyl compo- 
nents of the obtained metrics are presented. 

1. For the singular case we obtain three 5x5 non-singular matrices from 
(36). Let us denote \i k — -\- -\- z^.The first metric is given by (40). 
The non-zero Christoffel components are given 

pi _ pi _ pi _ pi _ pi _ r2 — r2 — _ r2 — r2 

24 — ^ 15 ~ '- 34 ~ ^ 35 ~ ^ 24 15 ~ 24 — o 34 ~ ^35 

y X y y y 

2 2 

_ _^p3 _ •^p3 _ _*^p3 _ p3 _ _'^'^p4 _ ^■^p4 _ ^ ^ p4 _ ^ p4 
— ^ 15 ~ ^ 24 — ^ 34 ~ ^ 35 ~ ^ 15 ~ ^ 24 — o i 2 34 ~ ^35 

z z y y y x^ + z'^ y 

2 

-Zl,,- ^124-^134-^2 + ^2^35- ^ 

Ricci scalar is given by i? = ^2_^_]p_^_^2 and Weyl tensor components are 

_ 12^^^ + 2y^ _ lyz _ ^ {^±j/)z ... _ 1 {z^ + y^)x 

H/1212 - -g , H/1213 - -3-^' '^1214 - -3 p , H/1423 - -3 p 

1 Ti/ 1 Ti/ 1 z{x'^ + z"^) 1 xyz 

yyi215 — W/1224 = 3~p~' '^1223 — 3"p"' '^1225 — "3 p ) 1^2335 — 3"^" 

•^1234 - 3-^, 1^1234 - 3-p-, 14^1313 - -3 p , 1^2323 - -3 ^ 

w -w -^""'y w -'^i^^y^w -w - ^"^y^ w - 

1^^1523 — l'V'2324 — 3"^"' ''^''^1314 — " ^ , I'1'1315 — I'1'1324 — "3"^") ^^''1323 — "3"^ 

_ 1 {z^ + x^)y _ 1 xyz _ 1 y^z _ 1 x^z _ 1 {x^ + z^)x 

1^1'1325 — 77 P5 , ^1^1334 — ~77~r9~' '''^1335 — " 77T9"' '^''2334 — 77-79"' '^''2325 — " 77 • 

6 k^ 6 k^ 6 k^ 6 k^ 6 k^ 

The second metric is given by 
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5(2) 



The non-zero Christoffel symbols are given by 





1 











~y 


\ 







1 





—z 


X 












1 


y 












—z 


V 










I 


-y 


X 











/ 



(63) 



pi _ y-pi 

^ 15 — ^ ^ 24 

Z 



X 



^ 15 



y^i 

x 



i34-^il5 



z 



^p2 
^24 



^p5 
^i24 



xy 



x + y 

5 _ „,T^5 



.r4 

24 



.r2 

^ 25 



l/'p4 



^p2 

y'' 



xy^4 

z'' 



^p3 



y^l. 



2 , 2^25 ~ 2/^35 

2;^ + 



a^p3 _ 



Iz2 + 2x2 + 2y2 _ l^y _ 1(^2 + ^2)^ 
3 P ' ^^^^ ~ ~3"fc2 ' ^^^^ " ~3 P ' 



Ricci scalar R 
The Weyl tensor components are 

W^1212 = 
W^1223 = 
W^1313 = 

W^1324 
W^1423 = - 



12 



^2 + y2 _|_ 



1215 



1 Z^X 

3l^ 



1 _ 1 xyz _ 1 2^2/ _ 1 z'^x _ 1 (x^ + 1/^)2: 

--p-, 14/1224 - --p-, H/1225 - --p-, ^Vl234 - --p-, H/1235 - p 



1 2z^ + 2x^ 
3 

1 y^x 



y2 _ 1 1/(2:^ + y2) ^ la;^^ ^ la^^ 

— , H/1314 - g p , H/1315 - ~3-p-' '^1323 - "g-p 



3 A;2 
1 (2;^ + 

'3 P ' '"^^^^ 

1 xyz 



A;2 

W^2325 



, W^1325 

,W^i, 



ly^ 

"3 F 

1 zx"^ 
3"P" 



, W^1334 = 

, W^2323 



1 xyz ^ 1 (x^ + y'^)y 

■3 F ' 3 A;2 



1 2^2 + 2?/2 + 



3 



2334 



3"P"' ^^2335 



A;2 

1 (x^ + y'^)x 
■3 P 



2324 



1 x^y 
3"P" 





1 








z 


-y 


\ 







1 








X 












1 


—X 














—X 










V 


-2/ 


X 











/ 



5(3) 



The non-zero Christoffel symbols of the second kind are given as 



(64) 



^pl _ '^p2 



r 



-r 

2/ 



25 



— r 



34 



^p3 
~Z^'' 



11 



'y 



^ 25 



^ 34 



2 14 — ,.15 



1/ 



•^^ 25 



JUL 34 



^p5 



p5 _ _^'^r5 _ ^^-pb _ _ZX 
2 15 — 25 — ^ 34 — ^ 



Ricci scalar R 



+ ~ y ' y 

12 



2 I o^2 I o„,2 



, Weyl tensor components are 



1212 



A;2 



H/1213 = -3^' I^ISU = 3-p-, 1^1215 = g-p- 



*^1223 — 3"p"' ^^1224 — g p , 1^1225 — 1^1234 — 3"^"' '^ISSS — " g p 

_ 1 22:^ + 20:^ + ^/^ _ ly^x _ Ixyz _ Ixy 

H/1313 - -g p , H/1314 - -g-p-' ^1315 - -3"^-' "^1323 - "g-p 



1324 



10x2+^ _ ly^z _ lzy\.. _l{x^ + f/)y 

3 P ' ^^^^ ~ ~3"P"' ^^^^ ~ ~3l^' ^^^^ ~ 3 



1423 



1 x'^y 
3^ 



1523 



1 zx^ 
" 3"fc2" 

1 xyz 



2325 



W^2323 
5 W^2334 = 



1 2^2 + + 



3 

1 xyz 



A;2 



2324 



A;2 

1 (x^ + 2;^)a: 
3 P 



3 ^2 ' " ^•^•5^ 3 ^2 ' " 2335 g ^2 

2. We are repeating the same procedure as above for the curved case. The 
metric given (55) has Ricci scalar is 0,VFi2i2 = Christoffel symbols of the 
second kind non-zero components 



1 (x^ + y^)x 
'3' 



pi __p3 _ x{l-y^) 2 __p4 _ 1/(1-^') 

^ 11 ~ ^ 13 ~ ' ^ 22 ~ ^ 24 ~ 



pi _ _p4 _ 
^ 22 ~ ^ 23 ~ 



x{l — X^) 



U 



p4 p3 p2 

^ 14 — ^ 24 — ^ 12 — 



U 



x"^ — y"^. 



xy 
u 



pi _ r3 — r4 — 

12 — ^ 23 ~ ^ 13 ~ 



X y 
u 



p2 _ p3 _ 
in— -l 14 — 



y(i-y') 



5 11 



U 



/ 1 + ^ 



u 
xy 
spa 

\ -y 



xy 
u 



1 + — — -v/u — X 



Ricci scalar is R=0,VFi2i2 = ~u 
The Christoffel symbols are 










J 



(65) 



pi _ x{l-y'^) 1 _ x'^y 1 _ x{l-x'^) 

1 11 — ) 1 12 — ) 1 22 — 

u u u 



p2 _ 1/(1 -y J T-2 
) ^ 11 



u 



1 ^ 12 

2 



xy 
u 



2 _?/(!- x') 



22 



U 



X 



^13 — ^14 — — ) ^24 — "V^J ^13 — /— 

y y 



, ^14^ — X 



12 
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Ricci scalar is R=0, VF1212 = — -• The Christoffel symbols are 



(66) 
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